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The displacement fluctuations of mirrors in optomechan- 
ical devices, induced via thermal expansion by temperature 
fluctuations due either to thermodynamic fluctuations or to 
fluctuations in the photon absorption, can be made smaller 
than quantum fluctuations, at the low temperatures, high re- 
flectivities and high light powers needed to readout displace- 
ments at the standard quantum limit. The result is relevant 
for the design of quantum limited gravitational-wave detec- 
tors, both "interferometers" and "bars", and for experiments 
to study directly mechanical motion in the quantum regime. 

PACS : 04.80.Nn, 05.40.-a, 42.50.Lc 



I. INTRODUCTION 

In a recent paper Braginsky et al JjJ , henceforth BGV, 
considered the noise in interferometric gravitational-wave 
detectors due to thermoelastic fluctuations of the mirrors 
attached to the test masses of the interferometer. These 
thermoelastic fluctuations have contributions from two 
independent processes, both acting via the thermal ex- 
pansivity of the mirror substrate material. The first one 
is the thermodynamic fluctuations in temperature of the 
body of the mirror substrate (these, in the approxima- 
tion of small thermal expansion, are independent from 
the thermodynamic fluctuations in volume, which are re- 
sponsible for the well studied thermal or brownian noise 
P|). The second one is the photothermal temperature 
fluctuations due to the fact that the number of photons 
absorbed by the mirror fluctuate. 

BGV results for the thermodynamic noise, obtained 
for half-infinite mirrors, have been extended to the case 
of finite size mirrors |3| , with particular reference to the 
design of advanced interferometric gravitational-wave de- 
tectors, such as LIGOII ||] . In both cases the calculations 
are concerned with mirrors at room temperature, made of 
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materials well in use for mirrors substrates as fused silica 
and sapphire, with km long Fabry-Perot cavities, which 
are characterized by laser beam spots of size r$ ~ 1 cm 
and comparatively low finesse T ~ 100, and with char- 
acteristic frequencies / ~ 100 Hz, for the mechanical 
motion to be monitored optically. 

There are a number of situations, at variance with 
the above, which are of interest for optomechanical de- 
vices. In such situations one or both of the thermoelastic 
fluctuations effects may be of concern, when one would 
like to reach, in the measurements of small displacement, 
the so called Standard Quantum Limit (SQL) P,k|. Al- 
ready for LIGOII, BGV seems to discourage, in favour of 
fused silica, the use of sapphire, which on the other hand 
may be the material considered for cold mirrors in con- 
nection with advanced configurations of interferometric 
gravitational- wave detectors, under study as LCGT M, 
LIGOIII @ and EURO §. 

The BGV effects would be of concern for very sensitive 
displacement sensors based on high-finesse Fabry-Perot 
cavities, to be used in connection with bar detectors of 
gravitational waves as dual cavity transducers PL or to 
study the quantum effects of radiation pressure [p_0| 12[ . 
In both cases the cavities are much shorter (less than 
a few centimeter) than in a gravitational-wave interfer- 
ometer, the beam spots are smaller (ro — 10 -2 cm), fi- 
nesses much larger {T > 10 5 ) and temperatures as low 
as T < 1 K. It may appear from BGV results that the 
thermoelastic effects would generate particularly large ef- 
fects, inasmuch the volume involved in the fluctuation 
processes would be correspondingly smaller. 

For these reasons, it is of interest to explore what would 
be the behaviour of both thermoelastic effects in the low 
temperatures and small beam spot regimes, where some 
BGV assumptions break down. In particular, the heat 
diffusion length l t depends on the temperature and can 
become larger than the laser beam spot dimension r , so 
that the adiabatic approximation is no longer valid. 

In Section y we gi ye the essentials of the regime of 
phonons and heat propagation, which establishes at low 
enough temperatures, and we evaluate the thermoelastic 
noises with a simple calculation, in relation to the beam 
spot size and to the frequencies at which the optome- 
chanical device is most sensitive. 

In Sections [II and |V we give an exact calculation of 
both thermoelastic effects in the whole region of interest, 



that is for any value of the ratio lt/ro- The results, under 
the assumptions of low temperature regime of Section ||, 
would directly apply to actual mirrors for the quoted 
optomechanical devices. We also relate in a general way 
the photothermal noise to the displacement noise induced 
on the mirror by the quantum fluctuations of radiation 
pressure in the cavity. 

In Section |V] we discuss limitations and relevance of 
our approach in the design of the SQL optomechanical 
displacement devices. 



II. HEAT PROPAGATION AT LOW 
TEMPERATURES 

Let us assume the optomechanical device to work in 
some frequency range centered around a frequency / and 
let us discuss the photothermal effect. We revisit the 
calculation of BGV in the following way, so to use it to 
see the regime which sets up at low temperatures. 

The multilayers coating of the mirror absorbs a small 
fraction of the light power and this induces an inhomo- 
geneous increase of the temperature of the bulk. The 
absorbed power is a Poisson distributed random variable 
(the statistics of the absorbed photon will be discussed 
in more details in section IV), and these fluctuations lead 
to thermal fluctuations in the bulk of the mirror. They 
are consequently responsible for fluctuations of the posi- 
tion of the reflecting face of the mirror, via the thermal 
expansion of the mirror material. 

The r.m.s. displacement noise of the mirror end face 
Az = zaAT is found by evaluating the r.m.s. fluctua- 
tion in temperature AT, in a volume V of the mirror of 
thickness z, linear thermal expansion coefficient a(T) and 
specific thermal capacity C(T), as the absorbed photon 
flux n fluctuates, 
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where Uojq is the energy per photon, An = y/ii/f is the 
r.m.s. poissonian fluctuation of the number of photons 
absorbed over the time 1/ / (n is the average absorbed 
photon flux), and p is the density of the mirror material 
(axis z is taken normal to the plane face of the mirror). 
At room temperature and for large beam spots, BGV 
conditions apply: the phonon mean free path and relax- 
ation times are very small respectively in comparison to 
the mirror coating thickness (where the photons create 
the phonons in the absorption process), and in compari- 
son with the characteristic time 1//. The thermal diffu- 
sion length at frequency / is given by 
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where K is the thermal conductivity. For a frequency / 
around fOO Hz, l t is on one hand larger than the coating 



thickness z c , and on the other hand much smaller than 
the beam spot radius rg, 



z c < l t < r Q . 



(3) 



This is the basic BGV approximation, which gives that 
the volume involved in the fluctuating thermal expansion 
effect is the fraction of mirror substrate V ~ Utq and thus 
one has z ~ l t in Eq. (|l|). 

This argument reproduces the essential features of 
BGV spectral density S z [/] of photothermal displace- 
ment noise, as one may write around the frequency /, 
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where S a b s — fiuoWafo is the spectral power noise of the 
absorbed light, with W a bs = hu^ri the average absorbed 
light power. In fact we see that Eq. (Q) is the same 
final BGV relation (Eq. 8 of H), apart from a term with 
the Poisson ratio of the mirror material and numerical 
factors. 

The condition (|3j) may break down for small beam spot 
radius r$ or for low temperature T , as the thermal diffu- 
sion length gets longer, cither in the mirror substrate or 
in the mirror coating or in both. 

For mirrors substrates of crystalline materials, as 
specifically sapphire, for a frequency / ~ 1 kHz, the 
thermal length l ts in the substrate at low temperature 
gets of the order of 10 cm, to be compared with a room 
temperature value /t s (300 K) ~ 10~ 2 cm (see Table |). 
Then at low temperature we rather have 
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at all frequencies below some I kHz, both for mirrors 
of gravitational-wave interferometers (for which rg ~ 
1.5 cm), and for optomechanical sensors (for which ro < 
3 10 -2 cm). This value for l ts stays constant in the whole 
region T < 10 K, as crystalline materials follow Debye 
T 3 laws for a(T), C(T) and n(T), and thus their ratios 
are all independent of T. 

High reflection coatings are typically 40 layers one 
quarter wavelength thick of alternating amorphous ma- 
terials as Ti02 and Si02, with a total thickness z c ~ 
10 -3 cm for Nd-Yag laser light. For such a coating, 
the breakdown temperature for Eq. (0) is different for 
mirrors of large-scale interferometers and for mirrors of 
high- finesse cavities, because of the difference in r$. For 
LIGOII mirrors for instance, taking Si02 as the refer- 
ence material for the coating, the thermal length l tc in 
the coating is of the order of ro only at very low tem- 
perature, T < 1 K. Amorphous silica films would have 
ltc > 10~ 2 cm for T < 10 K, so that for high-finesse 
cavities we have lt c ^ tq in the whole low temperature 
region. 

Despite this difference, there are two features of rel- 
evance, which affect similarly the thermal behaviour of 
the coating-substrate composite in both types of mirrors. 



Fused silica 


300 K 


10 K 


1 K 


a {K- 1 ) 


5.5 IO" 7 


-2.6 IO -7 


-2.6 10- 10 


K (W/m.K) 


1.4 


0.1 


2 10~ 2 


C (J/Kg.K) 


6.7 10 2 


3 


3 10- 3 


A (m) 


8 10- 10 


8 IO -8 


9 10- 6 


oi/k (m/W) 


3.9 IO -7 


2.6 10" 6 


1.3 10" 8 


h (m) 


3 10- 5 


1.2 1Q- 4 


1.7 1Q- 3 



Sapphire 


300 K 


10 if 


1 K 


a (K- 1 ) 


5 10- 6 


5.8 IO" 10 


5.8 IO" 13 


K (W/m.K) 


40 


4.3 10 3 


4.3 


C (J/Kg.K) 


7.9 10 2 


8.9 IO" 2 


8.9 IO" 5 


A (ro) 


5 IO -9 


2.2 IO" 3 


a/n (m/W) 


1.2 IO -7 


1.4 IO -13 


It (m) 


1.1 io- 4 




0.11 



TABLE I. Thermal properties of fused silica (top) and sap- 
phire (bottom) at different temperatures. The thermal expan- 
sion coefficient a, thermal conductivity k, thermal capacity 
C and phonon mean free path A are derived from |l| at room 
temperature, and from [[Cfl-llq] at low temperatures. The ther- 
mal length l t at 1 kHz is obtained from Eq. (@). 



In both cases we have that, at all low temperatures, the 
thermal length l tc stays longer than the coating thick- 
ness, Itc > z c , and that the mean free path A s of the 
phonons in the substrate is itself long, at least a fraction 
of cm |13| . In a coating of a Si02 film even the phonon 
mean free path A c will be larger than IO -3 cm, and thus 
\ c > z c ,forT<lK @. 

Let us then consider how the thermal regime changes 
at sufficiently low temperatures (T < 10 A"). The heat 
delivered by the absorbed photons in the volume r^z c 
of the coating crosses to the substrate in a time smaller 
than 1//, as lt c > z c . From the substrate, as A s > ro, the 
thermal phonons thereby created will reenter the coating, 
heating it up in even shorter time over distances A s . This 
happens because the acoustic mismatch between coating 
and substrate is small U%, when densities and sound 
velocities are quite close. The substrate thus acts as a 
thermal short for the coating in the plane of the mirror 
end face: the coating and the substrate will be practically 
isothermal over distances of the order of the phonon mean 
free path X s in the substrate. Then the coating will con- 
tribute to the thermoelastic fluctuations with its thermal 
expansion coefficient a c , but following the thermal fluc- 
tuations of spectral density St of the substrate. On the 
other hand, at the frequency /, the volume of substrate 
involved in the fluctuating heating will be of the order 
V — lt s , where the thermal length is that in the sub- 
strate. So, including both the coating and the substrate, 
we write now for the displacement spectral density S z [/] : 



S: 



[/] ~ ((a c z c ) 2 + (a s l ts ) 2 ) S T [f] 
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According to table |, a c z c is at least one order of magni- 
tude smaller than a s lts at low temperature (z c ~ IO -5 m 
and lt a — 0.1 ro). We can then neglect the expansion of 
the coating over its thickness z c and we find that the 
effect is dominated by the substrate properties, 
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This is the relevant result of our discussion of thermal be- 
haviour of the coating-substrate composite at low tem- 
perature, in that now the temperature fluctuations in- 
volve comparatively large substrate volumes, instead of 
the comparatively small coating volume, where the ac- 
tual absorption of photons occurs. Notice that, would 
not this be the case, one would have of course very large 
effects just concentrated in the volume, the external sur- 
face of which is that where displacements are going to be 
measured at SQL sensitivities. 

When we substitute in Eq. ([?]) the expression (|2|) for 
the thermal length, we see a dramatic change of regime, 
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where now the (substrate) thermal conductivity k s ap- 
pears, instead of the thermal capacity, and the frequency 
dependence has disappeared. In fact the system behaves 
as in the low frequency region of a low pass filter, while, 
under BGV conditions, it was rather in the high fre- 
quency region. 

We develop in the following sections a rigorous cal- 
culation of the effects in the low temperature regime, 
which gives in clear details the features grossly antici- 
pated above and which can be directly applied to mirrors 
of interest for optomechanical devices, when the above 
thermal behaviour of the coating-substrate system is re- 
alized. 



III. THERMODYNAMIC NOISE 

In this section we determine the thermodynamic noise 
without any assumption on the ratio lt/ro between the 
thermal diffusion length in the substrate and the beam 
spot size. Our analysis is an extension of the procedure 
developed by Liu and Thorne []3| , but it is valid even when 
the adiabatic approximation is not satisfied. According 
to Eq. (|2|), the condition l t < ro can actually be written 
as a condition over the frequency, / > k/ ' pCrf y Our 
treatment is thus also valid for an angular frequency w = 
2nf smaller than the adiabatic limit lu c defined as 



pCri 
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As shown in the previous section we neglect any ef- 
fects of the coating, taking only into account the thermal 



properties of the substrate of the mirror. We also neglect 
any finite-size effects since we have shown that the vol- 
ume of substrate involved in the fluctuating heating is 
usually smaller than the size of the mirror, even at low 
temperature. We thus approximate the mirror as a half 
space, the coated plane face corresponding to the plane 
z = in cylindrical coordinates. 

The analysis is based on a general formulation of the 
fluctuation-dissipation theorem, used by Levin Q to 
compute the usual thermal noise (brownian motion) of 
the mirrors in a gravitational-wave interferometer. We 
know that in an interferometer or in a high finesse Fabry- 
Perot cavity, the light is sensitive to the normal displace- 
ment u z (z — 0, r, t) of the coated plane face of the mirror, 
spatially averaged over the beam profile. This averaged 
displacement u is defined as 



u(t) = 



d r u, 



(z = 0,r,t) 
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To compute the spectral density S~ [u>] of the displace- 
ment u at a given angular frequency uj, we determine the 
mechanical response of the mirror to a sinusoidally oscil- 
lating pressure. More precisely, we examine the effect of 
a pressure P (r, t) applied at every point r of the coated 
face of the mirror with the same spatial profile as the 
optical beam, 
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P(r,i) = -^e- r >°cos(w£), 



(11) 



where Fq is a constant force amplitude. We can com- 
pute the energy Wdi SS dissipated by the mirror in re- 
sponse to this force, averaged over a period 2it /u> of the 
pressure oscillations. The fluctuation-dissipation theorem 
then states that the spectral density of the displacement 
noise is given by 
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where k B is the Boltzman's constant. This approach has 
been used by Levin to compute the brownian noise p8| . 
We are interested here in the thermodynamic noise so 
that Wdiss corresponds to the energy dissipated by ther- 
moelastic heat flow. 

The rate of thcrmoelastic dissipation is given by the 
following expression (first term of Eq. (35.1)ofRef. [pi): 
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where the integral is on the entire volume of the mir- 
ror and the brackets (...) stand for an average over the 
oscillation period 2ir/u>. 5T is the temperature pertur- 
bation around the unperturbed value T, induced by the 
oscillating pressure. Wdiss 1S then related to the time 
derivative dS/dt of the mirror's entropy, which depends 
on the temperature gradient V<5T. 



To calculate the rate of energy dissipation Wdiss, it 
is necessary to solve a system of two coupled equations, 
the first one for the displacement u (r, t) at every point r 
inside the substrate, and the second one for the tem- 
perature perturbation 6T(r,t). As the time required 
for sound to travel across the mirror is usually smaller 
than the oscillation period 2tt/u>, we can use a quasistatic 
approximation and deduce the displacement u from the 
equation of static stress balance |19] , 

V (V.u) + (1 - 2a) V 2 u = -2a (1 + a) VST, (14) 

where a is the Poisson ratio of the substrate (a is the 
linear thermal expansion coefficient). The temperature 
perturbation ST evolves according to the thermal con- 
ductivity equation fll9[|, 
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where a 2 = k/ pC and E is the Young modulus of the 
substrate (p is the density, C is the specific thermal ca- 
pacity). 

The solutions of Eqs. (jl4|) and (|15|) must also fulfill 
boundary conditions. If we approximate the mirror as 
a half space, the temperature perturbation ST and the 
stress tensor cr^ must satisfy the following boundary con- 
ditions on the coated plane face of the mirror, 

a zz (z = 0,r,t) =P(r,t), (16a) 

<?zx (z = 0, r, t) = a zy (z = 0, r, t) = 0, (16b) 

d{ST) 



dz 



(z = 0,r,t) =0. 



(16c) 



The stress tensor aij is defined in presence of changes of 
temperature as (see Eq. (6.2) of pjl) 



T ^ = -Jl^2a) a5T ^ 
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where the strain tensor uu is equal to i ( |^ + -^r- 



"v 10 c h^ "" 2 \d^~ + a^T, 
We solve perturbatively this system of equations at the 
first order in a. We first solve the static stress-balance 
equation at the zeroth order in a, neglecting the tem- 
perature term in the right part of Eq. ( |l4| ) and in the 
expression ([17]) of the stress tensor. The solution u^ ^ of 
this equation is well known (paragraph 8 of fly]). We 
then solve the thermal conductivity equation ( ]l5[ ) using 
as a source term the solution u^ ^ and we obtain the tem- 
perature perturbation ST 1 * 1 ' in the first order in a. The 
calculation of u^ ) , ST^ and finally Wdiss is done in Ap- 
pendix A. Using the results of this appendix, we show 
that Sc- [to] is equal to 
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where the integral I is given by 
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with fc^ = fc 2 + A; 2 and k 2 = k 2 ± + k 2 . 

We can express <S~- [w] as a function of an integral J [0] 
which depends only on a dimensionless variable O equal 
to oj/lu c , where to c — a 2 /r 2 , corresponds to the adiabatic 
limit (see Eq. |). We get 
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where J[Oj is derived from the integral / by the trans- 
formation of variables u = k±ro and v = k z ro, 
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When w>w c (i.e. ^> 1), we can neglect u 2 +i> 2 with 
respect to O in the denominator of the integral. J [Cl] can 
then be calculated analytically and we obtain 



J[0> 1] = I/O 2 



(22) 



Using this result and the definition of O, we finally show 
that <S~ [oj] is equal to 



S- [w > w c ] 
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(23) 



This formula is identical to the expression (18) of Ref. 0] 
and to the expression (12) of BGV |Q. 




Frequency ra/co c 

FIG. 1. Frequency dependence of the thermodynamic 
noise. The frequency ui is normalized to the adiabatic limit 
Lu c . The dashed curve corresponds to the adiabatic approxi- 
mation. 

For small values of O, the integral J [f2] can be com- 
puted numerically. The result is shown in figure [fl where 
we have plotted J [CI] as a function of O in logarithmic 
scale, for O between 10~ 5 and 10 5 . This figure shows that 



O = 1 is a cut-off frequency. For O > 1 the curve has 
a slope equal to —2, whereas for < 1 the curve has a 
smaller slope of the order of —1/2. In this low frequency 
range, the noise is smaller than the one which would 
be obtained using the adiabatic approximation (dashed 
curve in figure |l|) . 



IV. PHOTOTHERMAL NOISE 

We now briefly examine the case of the photothermal 
noise which exhibits somewhat a similar frequency be- 
haviour as the thermodynamic noise. We use the same 
method as BGV [Q to calculate the spectral density 
iS~ [tv] due to this noise but we do not make any adia- 
batic approximation so that the calculation is valid also 
for frequencies smaller than the adiabatic limit uj c . We 
then obtain: 
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K[Ci] = 
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When u> S> to c the adiabatic approximation is valid and 
the result of BGV should be recovered. Indeed, when 
Cl ^> l,we can neglect u 2 + v 2 with respect to CI and 
calculate analytically K [Cl] which turns out to be equal 
to I/O 2 . Using the definition of O, we finally show that 
S~ [to] is equal to 
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This formula is identical to Eq. (8) of BGV Q. 

For low values of Cl, K [Cl] can be calculated numeri- 
cally. The result is shown in figure 0. As in the case of 
the thermodynamic effect (figure [l]), = 1 is a cut-off 
frequency: for > 1 the function has a slope equal to 
—2, whereas for O < 1 the function has a much smaller 
slope and is almost constant. 

This result is in perfect agreement with the simple 
derivation made in section ||. The spectral density (J24J) 
can actually be written as 



S^[Lo] = ^{l + af(^) 2 S abs K[Cl] 



(27) 



which is similar to Eq. (ph at low frequency where 
K [O] ~ 1, apart from a term with the Poisson ratio. The 
frequency dependence of the photothermal noise then 
corresponds to a low-pass filter, with a cut-off frequency 
equal to the adiabatic limit uj c . At low frequency, that 
is when the thermal diffusion length l t in the substrate 




Frequency co/(o c 

FIG. 2. Frequency dependence of the photothermal noise. 
The frequency w is normalized to the adiabatic limit u c . The 
dashed curve corresponds to the adiabatic approximation. 



becomes larger than the beam spot radius tq, one has a 
dramatic change of regime and the photothermal noise 
is much smaller than the one which would be obtained 
within the adiabatic approximation (dashed curve in fig- 
ure |). 

Another important point for the realization of optome- 
chanical sensors working at the quantum level is to com- 
pare the photothermal noise to the displacements in- 
duced by the quantum fluctuations of radiation pressure 
of light. Quantum effects can be made larger than the 
usual thermal noise by decreasing the temperature and 
by increasing the light power. This would not be conve- 
nient for photothermal noise since the effect is propor- 
tional to the light power. BGV results show that for 
sapphire at room temperature the photothermal noise is 
of the same order as the standard quantum limit in an 
interferometer. In a high-finesse cavity, hopefully, these 
two effects are related to quite different photon statistics. 
As a matter of fact, the displacement u induced by the 
radiation pressure of the intracavity field is related to the 
intracavity photon flux N, 

U [u] = Xeff M Prad [u] = 2kk X eff [u>] N [u] , (28) 

where 2Tik is the momentum exchange during a photon 
reflection (k is the wavevector of light), and Xeff ls an ef- 
fective susceptibility describing the mechanical response 
of the mirror to the radiation pressure P ra d [fL2|- The 
noise spectrum Sc [w] induced by radiation pressure is 
thus proportional to the spectral power noise S cav of the 
intracavity light, which for a resonant cavity is E2j: 
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where uj cav is the cavity bandwidth, T is the cavity fi- 
nesse and W cav = hui N is the average intracavity light 
power. At low frequency (uj < to cav ) the intracavity 
photon flux corresponds to a super-poissonian statistics, 
the noise power being larger than the poissonian spectral 
density by a factor 2 T jit M . 



On the other hand, the absorbed photons always corre- 
sponds to a poissonian statistics, even if it is not the case 
for the intracavity photons. The spectral power noise 
Sabs of the absorbed light is given by 



S abs = TlUJoWa 
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where A is the absorption coefficient of the mirror (the 
average flux of absorbed photons is n — AN) . This effect 
cannot be understood within the framework of a corpus- 
cular model in which the photon absorption is described 
as a poissonian process: due to the super-poissonian 
statistics of the intracavity photons, one would find 
a super-poissonian statistics for the absorbed photons. 
One has to take into account the interferences between 
the intracavity field and the vacuum fluctuations asso- 
ciated with the mirror losses. This can be done by us- 
ing a simple model where the absorption is described as 
a small transmission of the mirror and where the ab- 
sorbed photons are identified to the photons transmitted 
by the mirror. One thus has a high-finesse cavity with 
two input-output ports and it is well known that the pho- 
ton statistics of the light either reflected or transmitted 
by such a cavity are always poissonian, for coherent or 
vacuum incoming beams |21| . 

Eqs. (E9[) and (pG) clearly show that both the radia- 
tion pressure effect and the photothermal noise are pro- 
portional to the intracavity light power W cav ; however 
the displacements induced by radiation pressure have an 
extra dependence on the cavity finesse T . The photother- 
mal noise can thus become negligible as compared to 
quantum effects for a high-finesse cavity. 

To perform a quantitative comparison between the two 
effects, we calculate the susceptibility Xeff defined by 
Eq. (Eq). We determine here the mechanical response 
associated with the internal degrees of freedom of the 
mirror, which are of interest for displacement sensors. 
We thus ignore the radiation pressure effects associated 
with the global motion of suspended mirrors, which are 
the dominant contribution to SQL at low frequency in 
gravitational-wave interferometers. We calculate the av- 
erage displacement u induced by the radiation pressure 
P ra d, assuming the mirror is a half space (z > 0). The 
normal displacement u z (z = 0,r,£) of the coated face of 
the mirror can be deduced from the results of paragraph 
8 of 111, 
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r — r 



W ^ 



2hkN(t) . ,>/",, 

«*(* = 0- r ' *) = "^T (1 - a*) J dr 

Using the definition (nfl) of u, we obtain 

The integral can easily be calculated by using a new set 
of variables u = r r' and v = r + r' . We finally get 



I -a 2 



Xeff M 



2nEr, 



(33) 



and the noise spectrum S-- [ui] induced by radiation pres- 
sure fluctuations is equal to 






2-kEctq I 



(34) 



where c is the speed of light. 

For all the displacement sensors considered in this pa- 
per the characteristic angular frequency u is smaller than 
the cavity bandwidth u> cav ■ The noise spectrum S-~- [u>] is 
consequently independent of u> and equal to 



<S~ \U! -C U) r 




IT 



h,u W cl 



(35) 



This expression shows that the radiation pressure effect 
depends on the mechanical characteristics of the sub- 
strate (E and a) whereas the photothermal noise (Eq. 
|27| ) depends on the thermodynamic characteristics of the 
substrate via the ratio a/n. At low temperature, K [fl] is 
of the order of 1 and the ratio a/n is constant and equal 
to 1.4 10~ 13 m/W for sapphire (see Table g). E is equal 
to 4 10 11 J/m? and a is equal to 0.25 so that the ratio 
between the photothermal and radiation pressure noises 
is of the order of 



S pt/ S rad „ g g 1Q U^JL 
u ' u JT 



(36) 



For a 1 ppm absorption rate (A = 10~ 6 ), a beam spot 
size r of 10~ 4 m, and a cavity finesse T of 10 5 , the 
photothermal noise is more than 4 orders of magnitude 
smaller than the radiation pressure effects of internal de- 
grees of freedom of the mirror. The photothermal noise 
is thus negligible as compared to quantum effects in op- 
tomechanical sensors. 

Note that this is not the case in gravitational-wave in- 
terferometers where ro ~ 10 -2 m and T ~ 100. The pho- 
tothermal noise is then 2 orders of magnitude larger than 
the quantum noise of internal motion. However, the inter- 
ferometer is not expected to be sensitive to this quantum 
noise since for suspended mirrors it is overwhelmed by 
the quantum noise associated with external pendulum 
motion. 



V. DISCUSSION AND CONCLUSION 

We have shown that both thermoelastic and photother- 
mal noises have a frequency dependence which looks like 
a low-pass filter: below a cut-off frequency u C) these 
noises are much smaller than the noise which would 
be obtained according to the 1/cu 2 dependence at high- 
frequency. 

First lines in tables y gi ye the values of the cut-off 
frequency lo c /2tt for fused silica and sapphire, and for 
a beam spot size rp of 1 cm (first table) and 10~ 2 cm 



(second table). The results show that u c is increased 
when the temperature decreases (3 orders of magnitude 
for fused silica and 6 orders of magnitude for sapphire 
when the temperature is reduced from 300 K to 1 K). 
If we consider a typical frequency lu/2it of 100 Hz, the 
adiabatic approximation is never valid for sapphire at low 
temperature, whereas it is valid for fused silica only for 
large beam spot size. 



ro = 10 2 m 


Fused silica 


Sapphire 




300 A- 


1 K 


300^ 


1 K 


uj c /2tt (Hz) 


1.5 10" 3 


4.8 


2 10 -2 


1.9 10 4 


n 2 j{Q] 


1 


0.74 


0.98 


1.3 10~ 4 


S^ (m 2 /Hz) 


2.7 10~ 42 


3.4 10" 45 


1.5 10 -39 


2.6 10" 49 


iS~ (adiabatic) 


2.7 10~ 42 


4.6 10" 45 


1.5 10 -39 


2 10- 45 



'■() 



10" 



w c /2tt {Hz) 

n 2 J{Q] 

S~ <m 2 /Hz) 

S~- (adiabatic) 



Fused silica 
300 K IK 



15 

0.51 

1.4 10~ 36 

2.7 10 -36 



4.8 10 4 

3.5 10~ 5 

1.6 10~ 43 
4.6 10~ 39 



Sapphire 
300 K IK 



2 10 2 

6.4 10" 2 

1 10~ 34 

1.6 10~ 33 



1.9 10 8 

1.4 lO" 10 

2.8 10 -49 

2 10" 39 



TABLE II. Results for fused silica and sapphire at differ- 
ent temperatures, for a frequency ui/2-k = 100 Hz and for a 
beam spot size ro = 1 cm (top) and 10 -2 cm (bottom). The 
thermodynamic noise S~- (3 r lines) is reduced as compared 
to its value obtained within the adiabatic approximation (i th 
lines) by a factor D. 2 J[Q] (2 nd lines). 

We first focus on the thermodynamic noise whose val- 
ues calculated from Eq. (pet) are shown in the third lines 
of tables ||. The noise is smaller than the one which 
would be obtained within the adiabatic approximation 
(last lines in tables O). The reduction factor, equal to 
1/£7 2 J[J7], can be as large as 10 4 for sapphire at low 
temperature with ro = 1 cm, and as large as 10 10 for 
ro = 10~ 2 cm (second lines in tables |J). 

We immediately see the impact for gravitational-wave 
interferometers (ro = 1 cm): for sapphire at low tem- 
perature, the thermodynamic noise is more than 4 or- 
der of magnitude smaller than for fused silica, so that 
the choice of the material at low temperature would be 
just the opposite than, as in BGV, at room tempera- 
ture. Furthermore, the thermodynamic noise at 100 Hz 
would be equal to 2.6 10 -49 m 2 /Hz for sapphire at low 
temperature, well below the noise at room temperature 
(2.7 10 -42 m 2 /Hz for fused silica). It is also well be- 
low the SQL limit due to the external pendulum motion, 
equal to 3.6 10 -41 m 2 /Hz for a mirror mass of 30 kg [nl. 

Similarly for optomechanical systems with smaller 
beam spot size (ro ^S 10~ 2 cm), the thermodynamic 
noise can be made as small as 2.8 10~ 49 m 2 /Hz by using 
sapphire at low temperature, to be compared to a noise 
larger than 10~ 36 m 2 /Hz at room temperature both for 



sapphire and fused silica. It is worth noticing that this 
very low value is partly due to the reduction factor asso- 
ciated with the non adiabaticity which is of the order of 
10 10 . This noise can be compared to the SQL limit due 
to the internal motion of the mirror, which is equal to 
fr\Xeff\ — 10~ 42 m 2 /Hz p2| . At low temperature, the 
thermodynamic noise is thus smaller than the SQL limit 
so that optomechanical sensors as in Refs. ||[L0| would 
be able to get to the SQL limit. 

Let us note that the thermodynamic noise for sap- 
phire at 1 K is mostly independent of the beam spot 
size ro: similar values are obtained for large spot sizes 
(2.6 10~ 49 m 2 /Hz for r$ — 1 cm) and small ones 
(2.8 1CT 49 m 2 /Hz for r = 1(T 2 cm). This result is 
due to the fact that, in contrast with fused silica, the 
adiabatic approximation is not valid for sapphire what- 
ever the beam spot size is, as long as it is smaller than a 
few centimeters. The non adiabatic condition u < lo c can 
actually be written as r < yj ' n/ ' pCto (see Eq. ||) . In this 
non adiabatic regime, we have shown that J [Q] evolves 
as 1/v" which is proportional to yfuT c and then to 1/fo. 
The thermodynamic noise is proportional to r^J [51] (Eq. 
|20| ) and is then independent of r^ . 

Similar results can be obtained for the photothermal 
noise. We have shown that as long as the adiabatic con- 
dition is not satisfied, the noise mostly depends on the 
ratio a/n (Eqs. || or g7(). In particular it does not de- 
pend on the frequency u nor on the beam spot size Tq. 
For sapphire at low temperature and for an average ab- 
sorbed power W a bs — 1 W of Nd-Yag laser light, the 
photothermal noise is then of the order of 10 -45 m 2 / Hz 
both for interferometers and optomechanical sensors, well 
below the SQL limits of both the external and internal 
motions. 

Let us finally note that the results obtained above ap- 
ply in detail to an actual mirror system when the condi- 
tions described in section O are fulfilled. In particular the 
interplay of the various characteristic lengths (phonon 
mean free path and thermal lengths at the frequencies of 
interest, both in the substrate and in the coating, beam 
spot, coating thickness and mirror size) must in the end 
allow, in some temperature range, that the thermal prop- 
erties of the substrate dominate. 

This may be not easy to achieve and thus our analysis 
may correspond to a somewhat idealized situation. At 
the lowest temperatures T < 0.5 K, the phonon mean 
free path in the coating gets of the order of its thick- 
ness (as for an amorphous silica coating, see for example 
pq|), while the time constants for phonon local equilib- 
rium continue to stay smaller than 2ir/u> both for the 
coating and the substrate. At first sight, it may appear 
that the equalization of coating versus substrate temper- 
atures will be even more facilitated. However at some 
intermediate temperature below 10 K the phonon mean 
free path in substrates like sapphire may get so long to 
exceed the dimensions of the mirror. In this case a more 
ad hoc model has to be considered, in which one speci- 



fies the details of the thermal link of the mirror to the 
main heat sink. Similar care should be taken if the ther- 
mal length at the lowest frequency of interest is of the 
order of the mirror size. Possibly in both cases the effect 
would be smaller than predicted in this paper, because 
the characteristic times for thermal equilibrium in the 
mirror volume would get even shorter. In any case, as 
to quote just one instance, phonon mean free paths and 
thermal conductivities (and thus thermal lengths), have 
strong dependencies at low temperatures on the level of 
impurities, so that each experimental configuration may 
be a case per se. 

In conclusion our results may be possibly of interest 
in two respects. First because they let see promising 
the use of low temperatures both for gravitational-wave 
interferometers and for optomechanical devices. Second 
because they may be of help to study, in actual experi- 
mental configurations for SQL conditions, the role of the 
thermoelastic effects in respect to choices of substrate 
materials, working temperatures, cavity finesses, mirrors 
losses, beam spot sizes and laser powers. 
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APPENDIX A 

In this appendix, we calculate the spectral density 
fy- [ui] (Eq. pit) of the spatially averaged displacement 
u induced by the thermodynamic noise. We approxi- 
mate the mirror as an infinite half space (z > 0). At the 
zeroth order in the thermal expansion coefficient a, the 
solution u(°) of the quasistatic stress balance equation 
(RlJs given by a Green's tensor (Eqs. (8.13) and (8.18) 
of jl9(|). The pressure P(r,t) applied on the coated sur- 
face of the mirror has only a component along the normal 
axis z, and we obtain the following expression for the dis- 
placement expansion 0W = V.uC), 



9 (0) (r,i) = 



2 (1 + a) (1 - 2a) 

E 



F cos (u>t) x 



(27T) 2 



l^r -tiz+ifti+tjt/) 



(37) 



with fc i 



Z?2 _i_ JL2 



To calculate the dissipated energy Wdiss it is useful to 
analytically extend the pressure-induced expansion O^ ' 
for negative values of z in such a way that it is an even 
function of z. Its spatial Fourier transform 0(°) [k, t] is 
then equal to 



e (0) [k, t] = 



4(l+q)(l-2g) 
E 



F cos(ujt) — e 



k 2 



-Ik 2 



(38) 



Ky -T Ky . 



with k 2 = k 2 

In the same way, we analytically extend the temper- 
ature perturbation ST in the half space z < in such 
a way that ST 1 - 1 ' (r,t) is an even function of z. Using 
the Fourier transform of the thermal conductivity equa- 
tion (|15| ) and the expression (|3g) of O^ ', we find that 
(5TW [k, t] is equal to 



e Mt +c.c, 



5T { V [k, t]=A [k] 
where the function A [k] is given by 
2 (1 + a) aT iu>k± 



(39) 



A[k) 



pC k 2 {a 2 k 2 + iuj) 



F e-*^ r o. (40) 



We now determine the thermoelastic dissipation Wdiss ■ 
The integral over z in Eq. ( |13|) is limited to the volume 
of the mirror (infinite half space z > 0). Since ST^ (r,i) 
is an even function of z, Wdiss can be written as 



£( I d 3 r(v5T^ 



(41) M 



where the spatial integration is in the whole space. Us- 
ing the Bessel-Perseval relation, we can express the dissi- 
pated energy Wdiss as a function of the temporal average 



of \ST^ [k,t}\ which is equal to 2 \A [k]| 2 (Eq 

n f d 3 k 
2T 
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Using Eq. 



(2tt) 3 
j), we finally obtain 
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where the integral / is given by 
f _ , d 3 k a 2 k 2 



(2tt) 3 k 2 (a 4 fc 4 + lu 2 ) 
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(42) 



(43) 



(44) 



This expression allows to determine the spectral den- 
sity iS-~ [lo] of the displacement u from the fluctuation- 
dissipation theorem (Eq. O). One gets the result given 
in the text by Eq. @: 
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